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In this paper, we study the quantum decoherence induced by accumulation of electron tunnellings 
during the quantum measurement of a charge qubit. The charge qubit is a single electron confined 
in coupled quantum dots. The measurement of the qubit states is performed using a quantum point 
contact. A set of master equations for qubit states is derived within a non-equilibrium perturbation 
to the equilibrium reservoir due to the electron accumulation between the source and drain of the 
quantum point contact. The quantum decoherence of the qubit states arose from the electron 
accumulation during the measurement is explored in this framework, and several interesting results 
on charge qubit decoherence are obtained. 
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I. INTRODUCTION 

Quantum decoherence 0, 0, is mainly induced by 
the interaction of a microscopic system coupled with 
its environment. The system-environment coupling re- 
sults in a non-unitary evolution of the system, which de- 
stroys the purity of quantum states and leads to infor- 
mation loss toward the environment. These issues have 
attracted much attention in the study of quantum com- 
putation and quantum communication in recent years 
1^1 . Some investigations of quantum decoherence have 
been focused on the measurement induced decoherence 
[E El B S El DjJ and the dynamical controls of the 
decoherence |l2l 1 1 3| . For the implementation of realis- 
tic quantum information processors, these investigations 
become the most champion works in the field. In the 
solid-state quantum computer with charge qubits, quan- 
tum measurement of the qubit states can be performed 
by coupling the charge qubit with a sensitive electrom- 
eter such as quantum point contacts (QPCs) 0] or a 
single-electron transistor 0,0]. In this paper, we shall 
focus on the charge qubit measurement using a QPC, and 
study the non-equilibrium dynamical effect of the QPC 
on the decoherence of the qubit states, here the charge 
qubit is a single electron confined in coupled quantum 
dots (CQD). 

Literaturely, to study the quantum decoherence in- 
duced by quantum measurements, an equilibrium ap- 
proximation is applied to the reservoir state 0, IT?! . 
For example, in the study of a charge qubit measured 
by a QPC (Fig. QJ, the electronic reservoir, the source 
and the drain of the QPC, is assumed to be macroscopic 
enough in comparison with the CQD. Then the thermal 
equilibrium of the reservoir is kept continuously through 
rapid relaxation processes. This condition of a perfect 
heat bath causes electrons tunneling such that the extra 
electrons arriving at the drain will flow back rapidly into 
the source through a close loop of the transport circuit. 
No extra electron accumulates in the steady reservoir. 



Practically, however, the condition of the perfect heat 
bath is not guaranteed at mesoscopic scale of the QPC. 
The electron accumulation (EA) in the source and the 
drain of the QPC may destroy the equilibrium of reser- 
voir and influence the outcome of the electrometer. As a 
result, the effects of EA may become significant to the de- 
coherence of the measured qubit. In this paper, we shall 
treat approximately this intermediate non-equilibrium ef- 
fect as a perturbation to the equilibrium reservoir. The 
quantum decoherence of the charge qubit due to the EA 
is then studied within this framework. 

The paper is organized as follows: The model of the 
charge qubit measurement and the perturbation scheme 
of taking the EA effect into account are presented in Sec. 

II. A set of master equations for the reduced density op- 
erator of the qubit is also derived in this section. In Sec. 

III. , the quantum decoherence of the qubit in the thermal 
equilibrium limit is illuminated according to the resultant 
master equations. An abnormal dependence of the qubit 
decay mode to the temperature and bias voltage is ob- 
tained. The EA effect on the quantum decoherence of 
the qubit states is explored in Sec. IV. The dynamics of 
the EA is analyzed, and the EA induced decoherence of 
the qubit is extensively discussed. Remarkably, we also 
found a decoherenceless EA effect in the low temperature 
regime. Finally, a summary is given in Sec. V. 



II. CHARGE QUBIT MEASUREMENT USING 
QPC 



The model we consider here is a single electron confined 
in a CQD (as a charge qubit) that is measured by the 
QPC detector, see Fig.^,. The Hamiltonian of the whole 
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system is given by HE SHE! 

H = Hs + Hb+H', (1) 
H s = E L c+c L +E R c+c R + n (c+c R + c+c L ) , (2) 
H B = ^^eiafai + ^e r a+a r , (3) 

/ r 

H' = (n-5Qc^c R ) 1 ^2(a+a r + a+a l ) . (4) 

lr 

Here, Hg is the Hamiltonian of the charge qubit in which 
(cl) and (cr) are the electron creation (annihi- 
lation) operators of the electron sited in the two dots 
labelled by L and R as shown in Fig. la. They sat- 
isfy the condition c\cl + c r cr — 1. The parameter 
flo is the electron tunneling amplitude between the two 
dots. Hb denotes the Hamiltonian of the QPC reservoir, 
which is decomposed into two terms for the source and 
the drain, respectively, and a ; + (a;) and a+ (a r ) are the 
corresponding electron creation (annihilation) operators. 
Since the presence of an electron in a dot close proximity 
to the QPC will cause a variation in barrier of the QPC, 
n — > il — 5f2, the interaction Hamiltonian H' describes 
the coupling of the quantum dots and the QPC in close 
proximity to the dot R. The information of the qubit- 
states can then be read out through the QPC current p|- 

There are two channels for electron tunneling in the 
QPC. One is the forward tunneling processes from the 
source to the drain, and the other is the backward tun- 
neling processes from the drain to the source, see Fig.^i. 
In the high-bias regime, magnitude of the electron tun- 
neling rate for the forward tunneling processes is much 
larger than the one in the backward tunneling processes. 
The contribution of the backward tunneling processes 
could be ignored as treated in |5[. However, in a non- 
equilibrium description, both contributions of the for- 
ward tunneling processes and backward tunneling pro- 
cesses should be included for an arbitrary bias voltage. 
The Hilbert space is then spanned by the following basis 

{\cj>i) ® \D (L n ,R n )) , \cj>i) ® \D (L m ,R m ))} , (5) 

where \<j>i) is an electron eigen-state in the CQD (we 
define \4>q) = \g) and \4>\) — |e) as the ground and 
excited states in the eigen-state representation, respec- 
tively), \D(L n , R n )}= a+ • ■ • a+ n aj l ■ ■ ■ aj n \B ) is a state 
with n electrons accumulated in the drain through 
the forward tunneling processes, and \D(L m ,R m ^ = 
af ■ ■ ■ a ; + cif 1 ■ ■ ■ df m \Bq) a state with m electrons accu- 
mulated in the source through the backward tunneling 
processes. The state \B ) is the reservoir (equilibrium) 
vacuum, and the operator denotes the annihilation 
of an electron from \Bq) [or the creation of a hole be- 
low the Fermi energy of the source (drain)], while the 
operator denotes the creation of an electron above 
the Fermi energy of the source (drain). Therefore, an 
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FIG. 1: (color online), (a) Schematic plot of a single elec- 
tron in a coupled quantum dots measured by the QPC. (b) 
Schematic diagram of the electron accumulation in reservoir, 
the first figure shows the forward tunneling processes (FTP) 
of electrons from the source to the drain which causes elec- 
trons accumulating in the drain, and induces a variation of the 
Fermi energy. The second one shows the backward tunneling 
processes (BTP) of electrons from the drain to the source and 
a corresponding variation of the Fermi energy. 

arbitrary state of the whole system can be expressed as 

\m) = X>>®(E E hin Rn {t)\D(L n ,R n )) 

+ E E K LmRm {t)\D{L m ,R m ))). (6) 

To study the charge qubit measured by the QPC, the 
electronic reservoir is usually assumed to be in the reser- 
voir vacuum state (a thermal equilibrium state). The 
reservoir vacuum state is uniquely determined for a given 
temperature and chemical potential. The reservoir time 
correlation function is usually calculated in this state. 
However, tunneling electrons from the forward tunneling 
processes will occupy the levels above the Fermi energy 
of the drain for a short period and disturbs the reservoir 
vacuum state, see Fig. This phenomenon, called the 
electron accumulation (EA) in QPC, should be taken into 
account, vice versa for the backward tunneling processes. 
Meanwhile, the EA may also be induced by the impurity 
of the reservoir or a circuit with a worse transport mo- 
bility (the imperfect condition of the reservoir). Litera- 
turely, the effect of the EA is ignored [1 17111 liH Hl|. 
One assumes that electrons created in drain (for the for- 
ward tunneling processes) will be forced into the circuit 
instantly, no electrons are accumulated in the drain tern- 
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porarily. Therefore, there is no any variation to the Fermi 
energy. The read-out current is completely contributed 
by the created electron of drain through the circuit, as 
a result in the macroscopic limit. However, for a meso- 
scopic reservoir, a variation of the Fermi energy can be 
generated by the EA. It is certainly interesting to study 
what the influence of the EA effect is on the decoherence 
of the measured qubit. 

In general, the EA should be considered by treating 
the reservoir as a fully non-equilibrium system. However, 
the reservoir is indeed an asymptotic stationary state at 
mesoscopic scale. This is because the external bias can 
be regarded as a bigger macroscopic reservoir compared 
with the QPC (see Fig. The external bias forces the 
QPC reservoir to evolve into an asymptotic stationary 
state in a short period. For such a case, we can treat ap- 
proximately the intermediate non-equilibrium effect as a 
perturbation to the equilibrium reservoir. Thus, the EA 
can be described simply by a variation of the chemical 
potentials in the source and the drain. In the forward 
tunneling processes, the tunneling of n electrons moved 
from the source to the drain results in an effective in- 
crease of the chemical potential in the drain p R — > /x^ = 
p R + SfiR (n, ft) and an effective decrease of the chemical 
potential in the source ph — > p~[ = //£ — SpL (n,ft), see 
Fig. Here, ft — l/fc^T is the inverse temperature. In 
thermal equilibrium state, the electron distributions obey 

for the source 

Obviously, at 



to be described by 



Fcrmi-Dirac function Fi = 
and F r — ' 



l+exp/3(e(— fiL) 

for the drain. 



l+oxp/3(e r -p H ) 

a given temperature, the chemical potential determines 
the number density of electrons N in reservoir. The re- 
lation of its inverse function can be solved. Furthermore, 
because the QPC is a 2-dimensional electron gas, the 
density of states, gL, for the source (ga for the drain) 
can be assumed to be energy independent Then 
the chemical potential can be expressed as pl,r{N , ft) — 
4 In [exp(/3iV / gL,R) — l] ■ The variation due to the EA is 
given by p R x(N±n/A, ft) - fi R>L (N , ft) = ±5p, R , L (n, ft), 
where A is the area of the QPC. One can find that up to 
the first order of n/A, 



Sp-R,L(n,ft) 



A 



9R,L 



(7) 



In addition, the bias effect on the EA should be taken 
into account. The external bias and the QPC-dots in- 
teraction determine electron tunnelings in the QPC. The 
current partially induced by the QPC-dots interaction 
records the qubit information. The bias effect simply 
forces the reservoir into a local equilibrium state. That is, 
the fluctuation of the chemical potential arisen from the 
QPC-dots interaction will be suppressed by the bias ef- 
fect asymptotically, pl (t)—p R (t) — > Vd- The EA should 
decay to vanish (a electron-release process, ER for short) 
by the bias effect. This ER process is also associated with 
the imperfect condition of the reservoir. To simplify the 
problem, the stochastic EA process (also including the 
effects of the bias and QPC-dots interaction) is assumed 



5p>R,L(n,ft,t) = 8p, R>L (n,ft)e ^ cos(u d t). (8) 

The exponential decay factor in Eq. JSJ comes from the 
bias effect. The detailed expression of the ER decay rate 
Fr.l (n) depends on the imperfect condition of the QPC 
associated with the bias effect. In general, T r _l are func- 
tions of the tunneling electron number n, and are differ- 
ent for the forward tunneling processes and backward 
tunneling processes. 1/T r> l (n) denotes the release time 
of the ER. The processes repeat continually through the 
transport circuit which is given by the oscillation func- 
tion in Eq. ©, where Wd is the frequency of the repeated 
processes. It can be estimated that LUd ~ v t jl with I the 
length of the QPC in the tunneling direction and Vt the 
electron drift velocity in the transport circuit. Thus, Eq. 
© describes properly the stochastic nature of the tunnel- 
ing processes embedded in the ER process and the cyclic 
motion of the external transport circuit. Obviously, the 
magnitude of the chemical potential variation decays to 
zero as the reservoir approaches to the local equilibrium 
state, namely, pl (t)—p, R (t) = Vd- Similar to the forward 
tunneling processes, the tunneling of m electrons moved 
from the drain to the source for the backward tunnel- 
ing processes leads to the chemical potential decreasing 
PR — * Pr = PR — 8p R (m, ft, t) in the drain and increasing 
//£,—> = Pl + 8fj,L(m, ft, t) in the source. Up to the 
first order of 5fj, Rt L (n,ft,t), the perturbed Fermi-Dirac 
function is given by 



F±{n,ft,t)ttF r<l ±6n R>L (n,ft,t) 



dFrJ 

Our.l 



(9) 



As a result, even if the EA disturbs slightly the equi- 
librium of the reservoir, it may induce decoherence to 
the measured qubit. This effect can be explored from 
the master equation of the measured qubit which has the 
form as 



d Pit) - 1 f 
dt ih 



H s ,p{t)\-Rp{t), (10) 
where p(t) is the reduced density operator of the qubit 



(Tr cM [p tot {t)\ + Tr 3(n) [p tot {t)]) . (If) 



Ptot(t) in Eq. Hlff) denotes the total density 
operator of the whole system, and the condi- 
tional partial traces are defined by Tr D ( n ) [O] = 
E Ln , R AD(L n ,R n )\d\D(L n ,R n )) and Tr fl(m) [6] = 
^ Lmi5m ( J D( J L m ,i? m )|0|5(L m ,i? m )}. Here, TV flW 
(Tr^( m ) ) means to sum over all the allowed energy 
levels of the reservoir with the condition of n electrons 
tunneling the QPC barrier through the forward (back- 
ward) tunneling processes and then accumulating in the 
drain (source). R in Eq. (|I0|l denotes the dissipation 
operator, the detailed form of which can be found 
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in the Appendix. The qubit decoherence due to the 
electrical measurement is governed by the dissipation 
term —Rp{t) in Eq. (|10(1 . This term is composed 
of the reservoir spectrum functions, which are the 
Fourier transformation of the reservoir time correlation 
functions defined by Tr D („) (£ , (m)) [f + {t)f{t')p tot ] and 

Tr D(")(5<'"))[/W/ + (*')/0tot]- Tne reservoir fluctuation 
due to the QPC-dots interaction is encoded in the 
reservoir operators f + {t)f{t') and f{t)f + {t'), where 
ft = ^2 lr e lty6lr ~ eiS>t / h a+ai is the electron tunneling 
operator in the interaction picture of the reservoir. 
The tunneling processes of electrons relative to the EA 
have also been taken into account in these functions 
by treating the EA-fluctuated reservoir state as the 
perturbed Fermi-Dirac function through Eq. J2J). As a 
result, the decoherence of the measured qubit due to the 
EA effect is manifested in the master equation through 
the first order perturbation of the reservoir spectrum 
functions, which is associated with the perturbed term 
^PR,L{n, P, t) °f the perturbed Fermi-Dirac 

function in Eq. (j5J). A perturbation scheme taking the 
EA effect into account is then given as follows 



Ptot(t) = Po,tot(t) + £pi,tot{t) + £, 2 p2,tot{t) + 



(12) 



where £ = Ujp, is a perturbation parameter characteriz- 
ing the EA effect with U = ((5fj, L (n, 0) + 5Jx R {n,0)) jn 
and p, — (fiL + pr)/2. The zero order term po,tot{t) in 
Eq. i|12|l represents the usual states without the EA. The 
higher order terms describe the contribution of the EA. 
By taking partial trace to Eq. I|12JI , we have the perturba- 
tion scheme for the reduced density operator of the qubit 



p(t) = p (t)+£pi(t) + 



(13) 



The QPC, as a sensitive charge detector, may be re- 
quired a large electron transparency and a relatively 
strong coupling strength to maximize the detector sensi- 
tivity in the current experiments [^. ll8llT9l |. However, we 
shall focus in this paper on the tunneling transparency 
regime (the weak dot-QPC coupling). This is because, 
as Gurvitz has pointed out, the occupation probabil- 
ity of a single electron in the coupled dot (as a qubit) 
can be traced through the QPC readout current in this 
regime 5j. To take the weak dot-QPC coupling into ac- 
count, the second order cummulant expansion technique 
0- llOl. fill Il7j is used. Here, the Markovian approxima- 
tion has been introduced in the technique. It should be 
noted that since the electron tunneling in the QPC is 
determined by the external bias and the QPC-dots inter- 
action, it results in the qubit dynamics with three time 
scales: the qubit decoherence time, the ER decay time 
due to the electron release by the bias, and the correction 
time due to the QPC-dots interaction associated with 
the bias, temperature and the electron accumulation. In 
literature, the validity of the Markovian approximation 



could depend on the time scale of the correlation time 
[20|. When the correlation time is much smaller than 
other time scales (in the Markovian regime), the non- 
Markovian memory effect in the qubit dynamics can be 
ignored. The qubit decoherence under the environment 
effect through the QPC-dots interaction is simply a de- 
cay process. Taking the electron accumulation effect into 
account, an extra qubit relaxation may be induced but 
the enhancement should not be so significant under the 
perturbative treatment. Therefore, the Markovian ap- 
proximation [2C| is still applicable for the derivation of 
the master equations here. 

As a result, the master equations for the measured 
qubit in our perturbation scheme are obtained as follows 



e ■■ ^ t p (t) = -iL D p (t)-X[q,[[G^, Po (t) 



e ■■ j t Pi(t) = -iL D p 1 (t)-\[q,([[G^,p 1 (t) 



(14) 



+ ^(t)[[G< 1 \iVo(t)]])]. (15) 



Eq. I|14|) is nothing but the rate equation of the qubit for 
the reservoir in the thermal equilibrium state without 
the EA effect. The EA effect on the measured qubit 
is described by Eq. JEJ. In Eqs. (|14ll5jl . L D is the 
Liouvillian operator for the qubit, A = TrgLgR/h, the 



double-bracket commutator 



A,B 



AB 



(ab)\ 



and J-(t) represents the contribution of the bias effect to 
the EA 



- f .s , ,s e n +re « g R (1 + e 

F{t) = cos (uj d t) — ,r 

1 + r 



(16) 

Also, the operator q and for k = 0, 1 in Eqs. (|14I15|) 
are defined as 



q = P0-P1-P2, 



Pi. 



(17) 



i=0,l,2 



and the operators Po,i,2 H0| are given by 



Pn 



n- T li 



Sfl cos 9 



i\9) (sMe) (e|). 



<557 sin , , , , - <5f2sin# . . , , , . 
Pi = e) (g\ , P 2 = \g) e , (18) 



where 9 = cos 1 [(Er — El) /j], Er,l denote the en- 
ergy for a single electron state in right and left dot, and 



7 = y 40q + {Er — El) is the energy difference of the 
two eigenstates of the charge qubit. The coupled dots 
structure is shown by 9. The operators Pi. 2 are responsi- 
ble for the inelastic excitation and relaxation of the elec- 
tron in the dots coupling with the QPC. The coefficients 
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G± \ and G±\ in Eq. I|17(l are the reservoir spectrum 
functions without and with the EA 

G { Z=9 i0) (eV d ), GM=-figM (eV d ) , 



G 



(0) 

+.1 



(eV rf -7), G« 



r<(o) _ 

Lr_i_ o — 



-5 (0) 

,g(°) (eF d + 7 ) , G% = JigM (eV d + 7 ) (19) 



G 



(o) 



(jr +,i 



and 



3 (0) W 



G 



(i) 



P9 {1) (eVd - 7) , 

'Hi 

, (20) 



G 



(i) 



+ ,'! 



V d ->-V d 



1 - e - ^ ' 



d_ 

dp 







1 - e-P x 



(21) 



where is the bias voltage and e the electron charge. 
For convenience, e = 1 is used hereafter. Also, No(t) in 
Eq. I|15l) is the number density operator of the tunnel- 
ing electrons for the reservoir in the thermal equilibrium 
state 0- 

Meanwhile, the qubit information can be read out 
through the QPC current I = where N = 

Er=o" Tr fl(") \Ptot{t)] -ELo mTr D('") [Ptotit)]. An extra 
minus sign in the second term of N is added because elec- 
trons in the backward tunneling processes tunnel back- 
wardly from drain to source. In our perturbation scheme, 
I = ftNo + ^ft Ai + C(C 2 ) + - •■, and iVi(t) is the number 
density operator of the tunneling electrons for the reser- 
voir fluctuated by the EA. The QPC current is governed 
by 



C a : ^A> (t) = -iL D N (t) ~ ^{ 



—(0) \ -| 

G P o(t)q + H.c.U, (22) 

L D JVi(t)-A{[g,([[G(°),iVi(t) 
+ T{t) \\G (1 \Wo(t) 



Mo) ^(i) \ \ i 

G pi(t)+F{t)G N (t))q + H.c.)j, (23) 



—(k) 



where G = Ei=o,iA G +l ~ G -l)P* for fc = 0, 1. The 
operator Wo(t) in Eq. corresponding to the zero or- 
der perturbation of noise spectrum [lot [TJ is determined 
by the following master equation 



-W (t) = -iL D W (t) - \{ 

G^po(t) + 2G~ (0) N (t) 



H.c. 



}.(24) 



In the derivation of Eqs. (|15I23I) . the ER decay rates 
Tr^ have been assumed to be independent of the tun- 
neling electron number n for simplification. The detailed 
derivation of the above master equations can be found in 
the Appendix. 



As we can see, the EA effect is described by Eqs. 
(|15I23[) . Eqs. I|14I22I) are simply the rate equations for 
the thermal equilibrium state without the EA effect. If 
the bias current flows from the source to the drain such 
that electrons accumulate in the drain, the EA induces 
a negative effective bias voltage —SVea- Then, an ad- 



ditional current I\ = ^Tr 



m (t) 



is induced from the 



drain to the source. If the bias current flows from the 
drain to the source such that electrons accumulate in the 
source, a positive effective bias voltage SVea is induced 
and an additional current flows from the source to the 
drain. Eq. (|23|) shows a modification of the first order 
perturbation to the transport current in the QPC. Eqs. 
HI 4I24|I are the main result of the theory, which will be 
used to explore the quantum decoherence of the measured 
qubit in the rest of the paper. 



III. DEPHASING AND RELAXATION IN 
THERMAL EQUILIBRIUM LIMIT 

As we have discussed in the previous section, the EA 
dynamics can be considered as a perturbative modifica- 
tion of the thermal equilibrium dynamics for the reduced 
system. Without considering the EA effect, the electronic 
reservoir is treated as the thermal equilibrium state, in 
which the ER process is assumed to respond fast and ef- 
fectively. The measured qubit in this case can be studied 
by the zero-order perturbation master equation Eq. I|14f) . 
When the EA effect also becomes significant, the qubit 
undergoes a local non-equilibrium process initially, and 
then approaches to a stable state asymptotically due to 
the bias effect. In this section we will study the quan- 
tum decoherence of the qubit in the thermal equilibrium 
limit. The EA effect on the qubit decoherence will be 
explored in the next section. 

In order to show the intrinsic measurement effect on 
the decoherence of the qubit, we concentrate first on 
the qubit with symmetric CQD, namely, both dots have 
equal energy levels En = El and the qubit state is char- 
acterized by 9 = it/2. The qubit with asymmetric CQD 
will be studied later. 

There have been several efforts contributed on this is- 
sues using the secular approximation or discussed in 
the zero temperature limit [l_0j . Here, an analytic dis- 
cussion is presented. To discuss the dephasing and the 
relaxation of the qubit in the measurement processes, a 
set of matrix elements for the reduced density operator pg 
are introduced : p Q . r = po, ee - Po, gg , Pa a = Po.eg + Po,ge 
and ip ,p = Po,eg ~ Po,ge, where p 0<i j = (i\p \j), and 
\h3 = 9) (l*i 3 ' — e )) are the ground (excited) state of the 
qubit. Eq. i|14fl can then be expressed as 

P0,r (*) = -To jr (p ,r (*) - P0,r (oo)) , 



To,r = VdG+!ai PO.r (oo) 



r (0) 

rr(o) : 



(25) 
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and 

pQ.d {t) = -^PO,p (t) , PO.p {t) = ~frPO,d (t) ~ r o,rPO,p (t) , 

(26) 

where % — t^QlQr (S£l) /%, and 



i=l,2 \ 
i=l,2 \ 



(27) 



with fc = 0, 1 



According to Eq. I|25|l . the relaxation rate of the qubit 
Applying the previous definitions in Eqs. I|19I21[) . 
to 



is r , r 

it leads 



rjd ( Vd sinh (iV d - 7 sinh (3j) 
cosh /3Vd — cosh ^7 



P0,s 



-7% 



(28) 

A plot of r . r and po.r (°°) is shown in Fig. [21 The 
solid curves denote the symmetric CQD. The relaxation 
rate Tq^ is a positive-monotonic-increase function of the 
temperature and bias voltage. Increasing temperature 
and bias voltage will enhance qubit relaxation, as we ex- 
pected. The asymptotic matrix element po,r (00) denotes 
an asymptotic stable state of the qubit. It is a nega- 
tive increase-monotonic function of temperature and bias 
voltage with values in the range — 1 ~ 0. Because the 
external bias can be regarded as a macroscopic reservoir 
(an effective heat bath) with respect to the QPC, it is ex- 
pected that there exists an analogy between biased QPC 
and the heat bath y|. A limited case is po,r (00) — * as 
Vd » 7. The qubit tends to stay in a completely ran- 
dom mixed state in the high bias limit. This property 
is similar to the thermal randomness caused by the heat 
bath. According to po.r (00) in Eq. I|28|) . the ground state 
occupation of the qubit shows a similar dependence on 
the temperature effect and bias effect. Furthermore, it 
can be checked that To.,, — > 77^7 in the zero bias and zero 
temperature limit. The relaxation due to QPC can not 
be stopped by turning off bias. The QPC will exhaust 
energy and information of the qubit, as pointed out first 
in |f. 

The qubit dephasing is governed by the coupled differ- 
ential equations (|26|) . the solutions of which are a linear 
combination of two decay modes 



cie r °.p* + c 2 e r °-p t 



with dephasing rates r o 



(29) 

= 1^0, r i \^0,r 

x y' 1 — respectively. The mean de- 

phasing rate is Po, r /2. The coefficients ci 2 in Eq. I|29|) 
depend on the initial condition of the qubit. We then 
have the result 





FIG. 2: (color online), (a) The decay rate Tq t of the mea- 
sured qubit in the equilibrium state. The solid curves are 
for symmetric CQD (7 = 2f2o) and the dash curves for asym- 
metric CQD (7 = 10f2o). (b) The asymptotic reduced density 
matrix p ,r (°°) of the measured qubit. Two branches corre- 
sponds to the symmetric CQD (solid curves) and the asym- 
metric CQD with 7 = 10f2o(the dash curves), respectively. 



Po.p (00) -> 0, p ,d (00) -> 0. 



(30) 



The qubit completely dephasing asymptotically. In addi- 
tion, two decay modes show different decay behavior in 
the regime KT 0ir /2 > 7. A plot of r^ p versus Vd is shown 
in Fig- 01 The curves with non-smooth turning point de- 
note that Tq has a transition from a degenerate regime 
HTq^/2 < 7, in which both modes have the same de- 
phasing rate To >r /2, to a non-degenerate regime along Vd 
axis. In the degenerate regime, only the real part of 

is plotted. In Fig. the fast decay mode, rj p , shows 
that the dephasing rate increases as increasing tempera- 
ture and bias voltage. In the high temperature and high 
bias voltage limit, it approaches to the linear dependence 
as the usual expectation. However, for the slow decay 
mode, r^p, Fig. ^jp shows an upside-down dependence 
on the temperature and bias voltage. At a given bias volt- 
age, the dephasing rate of the slow decay mode decreases 
with temperature increasing. It also shows that the de- 
phasing rate decreases with the bias voltage increasing 
at a given temperature. The minimum dephasing rate 
appears in the large bias voltage and large temperature 
limit. This phenomenon is very different from the behav- 
ior of the fast decay mode that we have known. Even for 
the case of high density of state (i.e. with small value of 
r y/hr]d), the phenomenon is still preserved. However, it 
should be noted that for a quantum information process, 
quantum computations must deal with all possible ini- 
tial conditions. Therefore, the decoherence of the qubit 
is characterized by the larger dephasing rate rj p rather 
than the slow one, r^~ p , unless an ancillary operation on 
the qubit can be added to limit the initial conditions to 
a sub-domain belonging to the slow decay mode. 



Next we shall study the decoherence of the qubit with 
asymmetric CQD. In this case, Eq. (11411 can be expressed 
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FIG. 3: (color online), (a) The dephasing rate of the fast 
mode of the measured qubit in the equilibrium state, (b) The 
dephasing rate Tg p of the slow mode for the measured qubit. 
The figures are plotted with the fixed parameter 'y/hrjd — Qq. 



Pm(*) = TPo, P {t) -VdCos9 casSCrp po,d(t) 

+ sin 9(G^ b -G^ aPo ,r (*))], (31) 

Po,r(t) = rj d sm9 cos 9G^ p ,d(t) 

+ sin 9 (G%-G^ aPo At))}, (32) 



PO.pit) = -^-A),d(*) - Vd 



Gi 0) cos* 



+ Gf„sm9 2 ]p , P (t)A33) 



where G^ 0) = G$J + G^q = F d coth (/9F d /2). This is a 
set of coupled differential equations. There is no analytic 
form to define the decoherence rate like the case in the 
symmetric CQD. However, in the large bias limit 7 <C 
Vd we can find from the first two master equations the 
relation po,d(t) /f>o, r (t) — — cot 9. Then the relaxation 
rate can be obtained 



r n 



G 



(0) 



sin 9 l G 



(0) 

+ .n. 



(34) 



It can been shown that three decay rates from master 
equations I|31I33[1 correspond to the roots of the following 
equation 



= y 3 -2r°y 2 



T 



y-VdG\ 



(o) / TSinJ 



ft 
(35) 

The mean value of these roots is 2Tq r /3. r is also 
plotted in Fig. [SJi. All curves show a positive- monotonic- 
increase with the temperature and bias voltage. This 
bias and temperature dependence of the mean decay rate 
approaches to be linear in the high bias regime. In the 
low temperature limit, it becomes linear for almost the 



whole range of the bias voltage. In addition, it can be 
found in Fig. |2K that the asymmetric CQD has a longer 
decay time than the symmetric CQD. It indicates that 
the influence of the device structure on qubit decoherence 
can also become significant in the low temperature and 
low bias regime. 

The asymptotic behavior of Eqs. I|31I33[) can be solved 
easily 



P04 M 

Po,r (°°) 



Po,p (°°) °> 
-7 (cosh/3V d 



cosh f3j) 



Vd sinh l3Vd — 7 sinh /3j 



(36) 
(37) 



Pq r (00) implicitly depends on the device structure 
through 7. The plot of p® r (00) in Fig. |2b shows that 
there are two branches corresponding to the symmetric 
(solid curves) and the asymmetric (dash curves) CQD, re- 
spectively. At a given bias voltage, the symmetric CQD 
is asymptotically forced into a more random mixed state 
in comparison with the asymmetric CQD. This is be- 
cause the device structure effectively modifies the bias- 
dependence of the qubit coherence. 



IV. EA EFFECT ON QUBIT 

Now, we turn to discuss the properties of the EA and 
the decoherence of the charge qubit due to the EA effect. 



A. Properties of the EA 



The EA can be characterized by the EA induced cur- 
rent ii (t) = j|Tr £Aq (t) . h (t) is plotted in Fig. H 

which is calculated according to the resulted master equa- 
tions. The result can be studied in detail by the anal- 
ysis of Eq. J33J). In the right-hand side of Eq. 
because the traces of —iLoNi (t) and commutators van- 



ish, the terms XT (t) ( G 
A G p^q + H.c. 



-(i) 



N (t) q + H.c. = Ki and 



= Kq completely determine the 

EA current. Also, the factor T (t) of K\ is exponential 
decay in time, and P \ (to) <C No {to). Thus, K\ is ini- 
tially dominant, and then Ko becomes significant as K\ 
exponentially decays to zero. The linear increase of the 
EA current in the beginning is governed mainly by K\, 
see Fig. 0] 

Analytically, No (t) can be solved in the asymptotic 
limit 



"-0,99 (*) = C g0 + c glt, «0,ee (*) = C e0 + C e it, 



Cel 



co(G^ 



G 



(0) 



G 



(o) 



(38) 
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r = n 

LfL 



stable 




linear increasing 



t (in units : 1/£2 Q ) 



30 



FIG. 4: The EA current h (t) in the QPC. The qubit is ini- 
tially in the \L) state without energy-level offset. The fol- 
lowing measurement parameters are used to plot figures : 

n = Q , so, = o.iOo, P = i/Oo , fj-L = foon , (Vd = ion ) 

and gL,ii = 2/fia. The initial input number of electrons in 
the source is set to be 100. There is no electron accumulation 
at initial time. Two cases with different ER decay rate are 
plotted. 



given in Eq. Q45[l. It is worth noting that the qubit den- 
sity pi >r (t) records the history of the EA effect. Even the 
accumulated electrons has been exhausted in this stage, 
pi jT (t) in Ko(t) performs as an EA background field and 
affects the electron tunneling in the QPC. Also, the time 
independent relaxation shift p\_ r (oo) leads to a constant 
transmission probability for the first-order electron tun- 
neling. Finally, the stable behavior of the EA current is 
reached, see Fig. 0] As a result, with the EA effect, the 
asymptotic current up to the first order perturbation is 



Z (oo) + Moo) 



\8Vl 2 



sny v d 



7 2 )sinh/3F d 



2 Vd sinh [3V d — 7 sinh /3-y 
+ 2\£V d (fl 2 - Q5Q (l + pl >r (00) cas6)) . 

(42) 

B. Decoherence induced by EA effect 



where c g o and c e o are simply constants, cq = 

— \£l 2 Vd/G^?\, the higher order contributions of 5VL to 
the coefficient Co have also been ignored. For the weak 
interaction coupling, it can be numerically checked that 
n o,ge (t) <C ^o,gg(ee) (t) m the asymptotic limit. We then 
obtain 



Tr 



-T 



(t) ((2\n 2 ) 2 v d t- 



constant J 



(39) 



In the beginning (the K± dominant regime), the ER de- 
cay is adiabatic, namely, T (t) « constant (= T'). It 
leads to 



h (t) « -.F'f (2\n 2 ) 2 V d t + constant. 



(40) 



This shows why the EA current is linear increasing in 
time, as one can see in Fig. 0] 

In the crossover regime in Fig. 0] the EA is slowed 
down by the bias. The maximum amount of accumulated 
electrons occurs at the valley of the /i-curves, in which 
the maximum reverse EA current is induced. Obviously, 
the smaller the ER decay rate T^l, the stronger the EA 
as shown in Fig. 0] 

Due to the bias, electrons accumulated in the QPC are 



completely exhausted through the ER process. Tr 



A 



becomes 



of Eq. 12: il) vanishes in this stage, and Tr 

dominant. The asymptotic current can then be solved 
analytically 



/i(oo) 



-Tr 



ZN^t) = £Tr[tf (oo) 



dt 

2\£v d (n 2 - nsn (1 + Pl . r (oo) cos 6)) , 



(41) 



where the higher order contribution from Sfl 2 has been 
ignored, and the first-order relaxation shift px, r (oo) is 



The qubit relaxation and dephasing due to the 
EA effect are studied in this subsection. Eq. i|15|) 
describes this decoherence process. The terms 



XF(t) G (1) ,iVo(t) 



,q 



K{ in Eq. JUJ is re- 



lated to K\ of Eq. (|23[l through the relation Nk 

12n=o n Pk^f( t ) ~ ELo^il'W' Tnis term arises f rom 
the EA, and the induced relaxation is suppressed by the 

bias. To understand this property, let us look at the 



asymptotic matrix elements (i K^^oo) ij first. Apply- 
ing the previous result Eq. (|38|) and the asymptotic anal- 
ysis, we obtain, up to the order 



9 |^i (oo) 



gO, 

'g £j(oo) 
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(43) 



where 



K g0 w xnsnr (t) V d 



(G?>) 



and G 



(0,1) 
+ ,(a,b) 



is defined in Eq. I|27|) . It can be checked 



that KgQ is positive because of (G+ \ — G+' a ) < and 



t(0) 



(i) 



Since p x (t ) < N (i ), 
before the bias effect becomes 



-1 < (G 
i(i\pi(t)\i) 

active (i.e. in the EA effect dominant regime). We then 
have 

(9 \pi {t)\g) « K go\r=i t + constant, 

(e \p\ (t) \ e) « — K b q\ :f _ 1 t + constant. (44) 

This linear time-dependence behavior indeed indicates 
an extra qubit relaxation due to the EA effect. The nu- 
merical results in Fig. [3^ (plotted by green and black 
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curves) roughly show this linear time-dependence of the 
relaxation. In Fig. [SJi, the solid (dash) lines correspond 
to the ground (excited) state, respectively. Instructively, 
a numerical result for Tl : r = is plotted in Fig. |3Jd. 
A perfect linear character is obtained. The induced re- 
laxation from the pure EA effect obeys the linear time 
dependence, as shown by Eq. (|44|l . 



pure EA effect 



Cl 













(a) 





50 



CL 



(b) 




Cl ( 

0.0 


A/v 


L. ° 


t = 


(d) 



t (in units : l/a c ) 



t (in units : ) 



FIG. 5: (color online). The time evolution of the reduced 
density matrix of the first order perturbation. The symmetric 
CQD is simulated. The qubit is initially in \L) state with 
£Tr[pi (t = 0)] = 0.01. The input initial electron number is 
100. The following measurement parameters are used : = 
Q , Stl = O.lOo, = , Pl = lOOfio, g L ,R = 2/fi , 

cjd = and pr — pl ~ Vd with Vd = and Yl,r = 0.15fio 
for red curve, Vd = 10f2o and Yl,r = 0.15f2o black curve, 
Vd = 10f2o and Tl,r = 0.3flo green curve, and Vd = Wilo and 
Tl.b — cyan curve, respectively, (a) The qubit relaxation. 
The time evolution of the matrix elements p\ g (pie) for the 
qubit in the ground (excited) state are plotted by solid (dash) 
curves, respectively, (b) The qubit relaxation under the pure 
EA effect (Tl,r = 0). (c) and (d) The qubit dephasing. The 
time evolution of the matrix element p\d and p\ v are plotted, 
respectively. 



According to Eq. 143|) , it is obvious that this relaxation 
process is slowed down by the bias because K g $ reduces 
to zero by the decay factor J-{t). This slow-down effect 
causes the qubit decay to the lowest energy state (in the 
first order perturbation) in the relaxation process, see 

Fig. Et- Then, the term A G {0 \pi(t) ,q = K' of 

Eq. H15|) becomes dominant. The qubit state undergoes 
a transition from the relaxation to excitation, see Fig.|3Ji. 
After the accumulated electrons exhausted, this excita- 
tion is completely governed by K' . In addition, the black 
and green curves in Fig. [SJi show that a stronger EA ef- 
fect (i.e. with smaller ER decay rate) forces the qubit to 
experience a larger relaxation shift (i.e. a lower energy 
state). The asymptotic relaxation shift pi tr (oo) can be 



easily solved 

Pl,r (°°) 



(e \p x (oo)| e) - (g \px (oo)\g) 
—7 (cosh/3Vd — cosh/37) 
Vd sinh pVd — 7 sinh f3j 



(45) 



which is independent of the ER decay rate. As a result, 
the total qubit relaxation up to the first order perturba- 
tion is 



Ptot.r (00) 



1 + 



-7 (cosh f3Vd — cosh ^7) 



V d smhpV d 
2 /! + < 



7 sinh /?7 



(A«L + Hr) 



1 



Ag R 



(46) 



The perturbed term comes from the EA effect. 

We may also discuss the low bias voltage limit, Vd <C 7. 
It can be checked that -> as V d -> 0. Thus, K[ 
in Eq. \ib\ vanishes, and the EA makes almost no effect 
on the qubit dynamics. The result is shown by the red 
curve in Fig.[SJi. The qubit only experiences a relaxation 
process without the excitation one. 

Next, we shall calculate the decoherence rate of the 
qubit. For a low bias voltage, because of — > as 
Vd — » 0, the dynamic structure of Eq. ijTil) and ((T3|) are 
almost the same. The decay modes of the qubit con- 
tributed by the zero order and the first order pertur- 
bation are closer. The EA effect can only enhance the 
qubit relaxation. Dynamically, the decoherence rate of 
the qubit does not be speeded up by the EA effect. Ac- 
cording to Eq. I|28|) . we obtain the total relaxation time 
Ttot,r and the dephasing time T to t,p of the qubit with the 
EA effect 



Ttnt.r — 



1 



r ,rl 



COsh/37-1 

. n n ) J-t0t,p — ^J-tot.r I*'; 

rjdjsmhp-f 



for the symmetric CQD in the low bias limit. 

With bias voltage increasing, the first-order qubit de- 
phasing rate can not be solved exactly. We study the 
qubit dephasing according to Eq. (|15fl , or in terms of the 
set of coupled differential equations for the symmetric 
CQD 



Pi,d(t) 
Pi,p (*) 
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Pl,P (*) > 

" ftPhd (*) ~ r ,rPl,}j (t) 



where 



Tl,r — 



(48) 



(49) 



Pl,d = Pi, eg + Pl,ge, tPl,p = Pi, eg ~ Pl,ge and IH 0:P = 

n o,ge— n-o.eg- The term Tg^pi p (t) in Eq. I|48l) denotes the 
qubit dephasing arisen from the QPC-dots interaction as- 
sociated with the bias, while the term Ti^ r J- (t) no iP (t) is 
resulted from the EA. Because ri jr is proportional to the 
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mean chemical potential p (= vl+vr ) ; T l r can be larger 
than r r for p > ^ 7. The qubit dephasing is mainly 
determined by the ER decay rate and the mean chemical 
potential, it does not so sensitively depend on the bias 
voltage. However, in the low bias limit, especially for 
V d < 7, G^i -> (i.e. ri, r -> 0). The qubit dephasing 
rate becomes much smaller. The numerical plot of the 
time evolution of the qubit off-diagonal matrix element 
shows the coincident result, as one can see in Fig- an d 
d. 




Decoherenceless EA effect in the low 
temperature regime 



FIG. 6: (color online). The EA induced relaxation rate, 
curves are plotted with different temperatures. 



The 



We have pointed out that the qubit undergoes an ex- 
tra relaxation under the EA effect. However, this EA 
induced relaxation can be suppressed in the low temper- 
ature regime. The symmetric CQD is used to check this 
result. The qubit relaxation rate in the first order per- 
turbation can be deduced from Eq. (|15[) 



ri, r (t) 



-f>l,r (t) 



(pl,r (t) ~ Pl,r (OO)) 



0,r 



/C(t)I\ r (50) 



where K, (t) = J- it) no, r (t) /(pi, r (00) — Pi, r (t)), and 
n , r (00) = rigG^/G^ 
(n =Tr 



has been used 



N (t ) ), and r 1>r is defined in Eq. gHJl. The 

relaxation rate Ti tT (t) is time dependent. Due to the EA 
effect, initially, the qubit severely relaxes with the relax- 
ation rate r 1>r (t) w r^ r /C (t) for p > V d . After the EA 
effect suppressed by the bias, the relaxation rate asymp- 
totically reduces to the zero-order r , r . Here, IC (t) is a 
fluctuant factor. The EA induced relaxation rate I\ r is 
plotted in Fig. EI It shows that Ti jr has a very different 
behavior from ro, r - For a large mean chemical potential 
(p 3> Vd) in which p is almost bias voltage independent, 
the electron source of the EA is excited from the en- 
ergy levels below the Fermi surface. The variation of the 
chemical potential SVea induced by the EA is limited 
by the chemical potential. The <5VEA-dependence of the 
qubit relaxation rate leads to a bounded phenomenon, as 
shown in Fig. |SJ 

The relaxation rate Tx >r also shows an anti-symmetric 
bias dependence (see Fig. |SJ. This feature can be under- 
stood as follows. There are two kind of electron-tunneling 
correlations, G+\ and G_ i} are involved. G^\ (G_\) 
corresponds to the electron tunneling that one electron 
is from the source to the drain at r (t) and the other is 
from the drain to the source at t (r). Explicitly, G^)_s i 
is associated with Tr *> \f? f T pi(t)] (Tr D n[f t f+ pi(t)]) 
for forward tunneling processes and Tr^ m [/t frpi(t)] 
(Tr/jm [ftfr Pi{t)\) f° r backward tunneling processes. 
The relaxation rate is related to G^\ and G_\ through 

the relation T k . r = rj d E;=i,2 ( G +l + / 2 with k = 



0, 1. For k = 0, the electron-tunneling correlation is sim- 
ply governed by the QPC-dots interaction and the bias, 
no EA involved. It can be checked that a relation of the 
bias symmetry holds 



GfAV d ) + G { %(V d ) Gf.{-V d ) + G^A-Vd) 



(51) 



The To,,, does not be changed by reversing the bias (Vd — ► 
—Vd), due to the fact that the qubit relaxation caused 
by the electron tunneling in the QPC only depends on 
the amplitude of the external bias voltage. However, the 
electron-tunneling correlation under the EA effect (k = 
1) obeys the relation of the bias anti-symmetry 

G%{V d ) + G%{V d ) G^j-Vd) + G^(-Vd) 
2 2 

(52) 

To understand this relation, let us look at the forward 
tunneling processes, in which an extra negative bias volt- 
age is induced by the EA, V d — > V d — SVea- According 
to the result in Sec. IV B, an effective relaxation of the 
qubit is induced by the EA effect. On the other hand, 
for a negative bias voltage, an extra negative bias volt- 
age induced by the EA leads to — V d — > — V d — SVea- 
The amplitude of bias voltage is increased. An effec- 
tive excitation of the qubit is induced by the EA effect. 
Accordingly, Eq. I|52l) indicates that the EA induced re- 
laxation rate of the qubit under the positive bias voltage 
(Vd) is equal to the EA induced excitation rate of the 
qubit under the negative bias voltage (—Vd). 

In addition, the relaxation process can be separated 
into two effective modes 



Ti, r = VdP( 



R(z + ) + R(z-) 



P(V d ±j) 



(53) 



where the function R(z) = coth(z) — z csch(z) 2 is a step- 
like function with bounded values ±i. The function pre- 
serves the bias anti-symmetry. As discussed in Sec. Ill, 
the relaxation depends on the structure of the CQD. It 
leads to two effective modes corresponding to R(z+) and 
i?(z_). In the low temperature regime, these two ef- 
fective modes separate. Eq. I|52|) shows that the qubit 
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excitation rate in the R(zJ) mode cancels the qubit re- 
laxation rate in the R{z + ) mode in the regime of Vd lim- 
ited by ±7 . Then a relaxationless EA effect occurs, 
as shown in Fig. However, in the high temperature 
regime, the CQD-structure dependence of the relaxation 
is suppressed, and two modes R(z±) become degenerate. 
No relaxationless EA effect occurs. As a result, Ti yr (t) 
reduces to Tq^ (i.e. ri >r — > 0) in low temperature regime. 
For the qubit dephasing, it can be checked that the de- 
phasing term Ti^F (t) n ^ (t) — > as T — » 0. The qubit 
is also EA-dephasingless in the low temperature regime. 



APPENDIX: DERIVATION OF MASTER 
EQUATIONS 



The derivation of the master equations I|14I15I22I23|I is 

presented in this appendix. The total density operator 
in the interaction picture of the reservoir satisfies 



dt 



Pi (*) 



1 

ih 



H S +H'i,Pl{t) , H' I = e 



n He m . 



In the weak coupling regime, the interaction Hamiltonian 
H' can be treated by using the second order cummu- 
lant expansion technique fK ITol ITTl Il7| . By taking the 
conditional trace Tr D (*0(.6C m )) shown in Sec. II and the 
second order cummulant expansion, the master equation 
becomes 



(t) = -i'Lnpf^ (t) - Rfvpffa (t) , (A.2) 



(n) 



V. SUMMARY 



We have developed a perturbative theory to study the 
EA effect on the decoherence of the charge qubit, a sin- 
gle electron confined in CQD, measured by the QPC. The 
contribution due to the EA is treated perturbatively. A 
set of master equations for the reduced density matrix of 
the qubit has been obtained in this perturbation scheme. 
By solving the resulted master equations, we obtain the 
decay modes of the dephasing and the relaxation of the 
qubit in the thermal equilibrium limit, and study the 
temperature- and bias-dependence of the dephasing and 
the relaxation rate. We find two kinds of decay modes for 
the dephasing in the electrical measurement processes: 
one decay rate increases as the increase of temperature 
and bias voltage, as expected; the other shows an ab- 
normal dependence. The qubit in the later decay mode 
preserves the longest decoherence time in the high tem- 
perature and high bias limit, its time scale is much longer 
than that of the usual mode. In addition, the EA proper- 
ties are studied extensively. The EA current is obtained 
according to the master equations. We find that the EA 
current varies linearly with time under the pure EA ef- 
fect, and is then slowed down into an asymptotic stable 
state by the bias. The qubit decoherence due to the EA 
effect is studied based on this analysis. We find that 
under the EA effect the qubit relaxation and dephasing 
rate are much larger than the ones in the thermal equi- 
librium limit, unless the bias voltage is turned off. Also, 
the master equations show that an extra qubit relaxation 
is induced by the EA effect initially, and then be sup- 
pressed by the bias. Asymptotically, the qubit will be 
forced into the state with a small relaxation shift which 
is independent of the ER decay rate. However, in the low 
bias limit, the qubit will not be affected by the EA. The 
qubit decoherence rate will not be speeded up by the EA 
effect. Finally, we find a dccohcrcncclcss EA effect in the 
low temperature regime. A bias anti-symmetry in EA 
processes suppresses the qubit decoherence. These deco- 
herent behaviors worth being explored in experiments. 



dt 



Ri 



f(b)P%) (*) = 72 



dr Tr 



£>(")(£><")) 



{ 



h j t 

Hi(t), [G (t, r) H'j{t)G (t, r)+ , pj (t)] ] } , 

(A.3) 



where p^ b) (t) 



T^DM(DM)[ptot(t)] and G(t,r) the 
propagator of the CQD. The decoherence of the measured 
qubit is governed by the dissipation term — Rf{b)P^f%\ if) 
in Eq. (Ojl . 

Note that the master equations (|A. 21 IA.3J1 essentially 
describe a Markovian process. The Markovian approx- 
imation has been introduced in this stage. Therefore, 
the time integration in Eq. (|A.3|I is replaced by the one 
be integrated out along the whole time domain. Carry- 
ing out the commuter in Eq. l|A.3j) . Rf(b)P^fn,) (*) can be 
expressed as 



1 

2T? 



2h z 



dr 



{q(h 



■ 6 P\ 



x (C}\ bh _(t-r) + C}\ bh+ (t-r))} + 



dr 



{(Po 



r-«) 



H.c. 



-i 7 (T-t) » 

Pi - e s P 2 ) 



-7, + (b-) 



f,-(b,+y 



-))q}+H.c, 
(A.4) 



where the reservoir time correlation functions are defined 
by 



'/.+ (* 



Tr 



n -nftfrPdt)], 
+ , 



c 

Cl_{t- T ) = Tr DM [fJ?Pi(t)\ 
for the forward tunneling processes, and 

CT + (t-T) = Tr Dlm) [f+f TPl (t)], 



cr(t 



Tr 



[ftf?Pl(t)l 



(A.5) 



(A.6) 
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for the backward tunneling processes. It can be easily 
checked from Eqs. I|A.5IA.6(I that without the EA effect, 
the backward tunneling processes is covered in the for- 
ward tunneling processes, and the result is the same as 
that in |5j- 

Next, the EA effect is taken into account to calcu- 
late the reservoir time correlation functions. It should be 
noted that the excitation rate of an electron occupying 
in lower energy levels far from the Fermi energy is much 
smaller than the one near the Fermi energy. That is, most 
likely, only the electrons occupying near the Fermi energy 
can tunnel through the QPC. If the QPC has a low trans- 
mission, i.e. n/A, m/A < N, the states {\D (L n ,R n ))} 
contributed by all allowed tunnelings in which electrons 
occupying near the Fermi energy can tunnel from the 
source to the drain are almost the same. Therefore, we 
have the following approximation 

Tr D( Mft + frPi(t)} « Pf E (D n \f7fr\D n ), (A.7) 

L n R" 

and CJ + (t — t) can be expressed as 
#°(*)£l»H»^[^#MA where A n = 
• • • a£ ai ■ • ■ &i , and p^' is the reservoir vac- 
uum state. Furthermore, we define p^ = A n p^ A n + as 
the electronic reservoir density operator with n electrons 
created in the drain for the forward tunneling processes. 
Combine these analysis together, we obtain 

Cl + {t-r) = pf»(t)EE ef( ^ V)e ^ ( ^ r) 

V r' lr 

x Tr \a/^ aip^hi/, a r > 

lr 

x [l-F+(n,f3,t-r)], (A.8) 
where the EA-fluctuated Fermi-Dirac functions F^ r have 



been treated perturbatively and shown in Eq. (|5J. Note 
that Eq. © is valid when n/A <C N and m/A < 
N, which coincides with the approximation used in 
Eq. I|A.7|) . Other reservoir time correlation functions can 
be similarly reduced to the following forms 

Cl_(t -t) = ^W^e-^^-^F+tn,/?,^ r) 

lr 

x [l-F ; -(n,/3,i-r)], 
C£+(* - r) = p ( r\t) E e^^F+in, 0,t-r) 

lr 

x [l-F-(n,(3,t-r)], 
C™-(t - r) = p { b m \t) E e-^to-^F-M, t r) 

lr 

x [1 -Ff{n,0,t-T)]. (A.9) 

Obviously, it shows that, in the above time corre- 
lation functions, the perturbed terms proportional to 
3pn,L(n, f3, t) come from the EA effect. Using the pertur- 
bation scheme in Sec. II, we then have pj(t) = + 
^Pi^(t) + £?p2 f(t) + ■ ■ ■ for the forward tunneling pro- 

cesses, and p^\t) = j$(t) +tf§\t) + £ 2 # (*) + ■ ■ • 
for the backward tunneling processes. 

In addition, the time integration in Eq. I|A.4|) can be 
calculated by using the results 

1 f°° (* ) 

— J dre :L -^Cy (b)t+ (t - t) = 'KgL9Rpfl ) {t) 

x (.9 (0) (K i - 7 )T^5 (1) (Ki-7)) , 

— J dre"* *■ T C] {h) _(t -t)= irg L gR^ b) (t) 

x (g { < ) \-V d + 1 )±nUg {1 \-V d + 1 )) ,(A.10) 

and the relation 5pL (n,^,t) + Spa (n, (3,t) — np^F{t), 
the master equations become 



'dt 

x{q(Af + A^)p$(t) (A^p^\t) + A^p^m 

+ f[$(A?> + AW)p$(t) (Afp[y\t) + A^p^\t))q 

+ Ht)m{A ( i ) +AV)p$(t) -^((n-yAVp^it) + (n + 1)A^ p^(t))q] } + H.c 

o(e) + --- (a.ii) 



for the forward tunneling processes, and 
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|/«2(*)+€|/4S(t) 



in). 



+ il o) )pS (*) - Wr 1} w + 4 0) & +1) (*))« 

- ^(t)ng(i« + A^/gCt) + ^(t)((n - l)!^^) + (« + l^M? 



l \t))q])+H.c. 

(A.12) 



for the backward tunneling processes, where A 1 ^' 1 ^ = 

J2i=o i 2 G± 'P Pi- All the definitions of the elements of 
Eqs. (|A.lllA.12|i c an be found in Sec. II. 

Eqs. I|A.11IA.12|I describes the transport properties of 
electrons in the system we concerned, for example, the 
transport current and noise spectrum 0) B H DU DjJ • 
According to Eq. (JSJ and recalling that pffi b \(t) de- 
scribes the quantum oscillation of the electron in the 
CQD with the conditions of n electrons accumulating 
in the drain (source), the fc-th order perturbation of 
the qubit oscillation is completely described by Pk(t) = 

E^opJ/W + Em=opff (*)• Accordin ^y« the fc ; th or - 

der perturbation of the QPC current operator is given by 
where N k = E^fiQ® ~ E~=o^ff(*), and 



the fc-th order perturbation of noise spectrum is defined 
by W k(t) = EZ o^ (pj$(<) + Pklit)) ■ Combine with 
Eqs. HA.lllA.12p . one can obtain the master equations 
shown in Sec. II straightforwardly. 
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